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ABSTRACT 
he  PERT/COST  system,  now  in  current  use  a]  step: 

ing  assumptions  in  order  to  permit  the  use  of  analytic  methods  for  deter- 
??.ining  predictions  of  critical  path  activities  amid  project  duration.   As 
a  result  of  these  assumptions,  the  predictions  is  ;-.  involve  significant 
error,  and  will  always  be  optimistic.   The  Monte  Carlo  technique  of  net- 
work calculation  does  not  require  these  assumptions  and,  foence,  is 
capable  of  yielding  snore  accurate  predictions  and  providing  more  useful 
information.   This  technique  is  discussed  in  detail,  k   more  flexible 
probability  density  made!  for  activity  tisr»es  is  introduced,  and  a  re- 
source allocation  technique  based  on  the  probability  that  an  activity 
will  be  on  the  critical  path  is  developed.   Finally,  an  application  of 
PERT  network  theory  to  military  operational  planning  is  described. 
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1.    Introduction. 

The  PERT  system  was  devised  by  a  study  group  under  the  direction  of 
the  Special  Projects  Office,  Bureau  of  Naval  Weapons,  Navy  Department, 
in  order  to  develop  a  methodology  for  providing  the  management  of  the 
Fleet  Ballistic  Missile  (Polaris)  program  with  an  information  reduction 
system  for  program  monitoring  and  evaluation.   By  this  system,  management 
was  to  be  continuously  apprised  of  progress  to  date  for  the  program  as  a 
whole,  and  be  furnished  valid  predictions  as  to  outlook  toward  accom- 
plishing program  objectives.  [2J 

The  system,  as  developed,  was  a  giant  stride  in  management  tech- 
nology, and  is  credited  with  a  major  contribution  toward  the  rather 
phenomenal  success  of  the  FBM  development  program   Not  only  did  the 
PERT  system  accomplish  the  objectives  already  mentioned,  but  it  allowed 
management  to  predict  those  activities  in  the  development  project  whose 
completion  time  would  have  a  direct  effect  on  the  overall  project  dura- 
tion.  These  critical  activities  could  then  receive  appropriate  mana- 
gerial attention. 

The  PERT  system  has  since  been  employed  extensively  in  industry, 
and  has  been  made  a  standard  procedure  for  monitoring  research  and 
development  projects  under  the  cognizance  of  the  Federal  government. 
Industry  acceptance  of  the  system  has  been,  in  general,  enthusiastic  and 
widespread. 

The  system  has  been  recently  extended  to  include  a  cost  control 
feature,  PERT/COST,  which  has  now  been  designated  the  standardized  sys- 
tem for  monitoring  federally  sponsored  R&D  projects. 

Another  system,  the  Critical  Path  Method,  has  paralleled  PERT  in 


development.   It  is  similar  in  many  respects  to  PERT,  lacking  PERT's 
stochastic  representation  of  time,  but  incorporating  a  correlation  be- 
tween time  and  cost,  enabling  management  to  schedule  optimally.   The 
Department  of  Defense  and  NASA  have  incorporated  the  desirable  features 
of  both  systems  into  PERT/COST. 

The  success  of  PERT  and  Pert-type  systems  has  stimulated  the  inter- 
est of  professionals  in  the  fields  of  mathematics,  management  science, 
and  operations  research.   The  professional  journals  have  featured  many 
papers  on  the  relative  merits  of  the  systems ,  proposals  for  refinement 
and  extension  of  the  systems,  consolidation  and  integration  of  the  sys- 
tems, and  the  like.   It  is  the  purpose  of  this  paper  to  consolidate  many 
of  these  ideas,  and  to  propose  some  refinements  and  extensions  to  the 
system,  designed  to  extend  its  usefulness  and  increase  the  accuracy  of 
predictions  about  the  projects,  and  to  generate  specific  recommendations 
for  control  of  time  and  cost. 

We  shall  first  briefly  describe  the  basic  PERT  system,  critically 
discuss  certain  features  of  the  system,  and  tender  proposals  by  which 
the  basic  system  may  be  improved.   In  particular,  we  shall  discuss  the 
Monte  Carlo  method  of  calculation  of  the  PERT  network,  and  show  how  this 
method  overcomes  many  of  the  deficiencies  inherent  in  the  analytic  method. 

Next  we  shall  describe  the  standard  PERT/COST  system  currently  being 
implemented  by  DOD  and  NASA.   We  shall  follow  this  with  our  proposal  for 
an  integrated  project  control,  scheduling,  and  resource  allocation  tech- 
nique. 

Finally,  we  shall  discuss  an  extension  of  the  PERT  network  approach 
to  the  problem  of  planning  projects  involving  major  future  decisions, 
such  as  military  campaigns. 


2.    Conclusions  and  Recommendations 

As  a  result  of  this  study,  it  is  concluded  that: 

(a)  The  assumptions  employed  in  the  currently  standard  analytic 
PERT  network  calculations  may  induce  significant  error  in  the  results. 

(b)  The  standard  analytic  technique  for  network  calculation  fails 
to  provide  a  sufficiently  definitive  measure  for  activity  criticality. 

(c)  A  more  flexible  model  for  activity  completion  time,  providing  a 
range  of  variance  to  allow  for  varying  uncertainty  in  time  estimates,  is 
needed. 

(d)  The  Monte  Carlo  technique  does  not  depend  on  the  assumptions 
which  induce  the  inaccuracy  in  the  analytic  technique,  hence  is  capable 
of  yielding  results  limited  in  accuracy  only  by  the  validity  of  the 
estimates  for  activity  time.   In  addition,  this  method  yields  the  proba- 
bility that  a  given  activity  will  be  on  the  critical  path,  and  other 
useful  information  not  obtained  by  the  analytic  technique.   The  Monte 
Carlo  technique  can  easily  provide  the  flexibility  discussed  in  (c), 
above . 

(e)  The  Monte  Carlo  technique  requires  increased  computer  running 
time,  but  the  increase  is  not  prohibitive,,  and  is  well  justified  by  the 
improvement  in  quality  of  the  results  obtained, 

(f)  A  resource  allocation  technique,  based  on  the  probability  that 
an  activity  is  on  the  critical  path,  is  feasible,,  and  provides  a  reason- 
able method  of  optimum  resource  allocation. 

(g)  The  network  representation  of  the  complex  interrelationships 
between  activities  in  a  project  is  a  technique  which  may  well  be  applied 
to  the  task  of  planning  military  operations. 


It  is  recommended  that  the  Monte  Carlo  technique  for  PERT  network 
calculations  be  employed  in  lieu  of  the  analytic  technique.   A  method  of 
utilizing  this  technique  is  discussed  in  detail  in  Appendix  I. 


3.    The  Basic  PERT  System 

The  basic  PERT  system  is  a  means  of  defining  the  relationship  be- 
tween the  various  activities  comprising  a  project »  and  for  estimating 
the  time  required  to  complete  each  activity  aired  the  project  as  a  whole* 
Use  of  the  system  enables  realistic  schedules  to  be  imposed  on  the 
project,  and  assists  in  controlling  the  execution  of  the  project. 

By  a  project  we  shall  mean  a  collection  of  tasks,,  or  activities , 
each  a  necessary  step  toward  the  achievement  of  some  final  objective. 
Subsets,  called  paths,  of  these  activities,  are  dependent  sequentially s 
that  is,  any  activity  in  the  path,  other  than  the  beginning  activity 9 
may  not  be  commenced  until  its  predecessor  has  been  completed.   Each 
activity  in  the  path,  except  the  beginning  activity,  has  a  predecessor 
in  the  path.   An  example  of  a  project  in  this  sense  may  be  as  simple  as 
the  construction  of  a  house,  where  the  activities  are  tasks  such  as  lay- 
ing foundation,  erecting  subfloor^,  installing  rough  plumbing,,  erecting 
framing,  installing  rough  wiring,  etc.   Or  a  project  may  be  as  complex 
as  development  of  the  Fleet  Ballistic  Missile  system 9  with  many  hundreds 
of  activities  involved. 

A  network  of  directed  arcs  is  a  very  convenient  way  of  representing 
a  project.   The  activities  are  represented  by  directed  line  segments 9 
called  arcs,  terminating  at  nodes,  called  events.   (See  Figure  1.)  The 
time  required  to  complete  an  activity  is  associated  with  the  length  of 
the  arc  representing  it.   Events  represent  points  in  time.   Several 
activities  may  terminate  at  one  event.   If  so,  the  event  is  said  to  occur 
when  the  last  of  these  terminating  activities  is  completed.   At  that 
point  in  time,  other  activities,  whose  commencement  is  contingent  upon 


Beginning 


Figure  1 
A  Network  Representation  of  a  Project 


the  event,  may  commence.   These  activities  are  represented  by  arcs  di- 
rected outward  from  the  event.   In  the  house  construction  example ,  the 
installation  of  the  interior  wallboard  may  commence  only  after  the  rouj 
plumbing  and  rough  wiring  are  completed. 


Install       interior  ^_ 

VV<3)  I  boa  fd 


Figure  2 
Two  Independent  Activities  Incident  to  an  Event 

A  project  can  always  be  represented  by  a  network  with  one  beginninj 
event  and  one  end  event.   Fictitious,  or  dummy ,  activities  may  be  drawn 


in  the  network  In  order  to  achieve  this  format.   Such  dummy  activities 
may  consume  zero  time,  or  may  represent  a  predetermined  amount  of  sched- 
uled "dead"  time.   It  is  coM»n  practice  for  the  network  to  be  drawn 
with  beginning  event  on  the  leftv  and  with  the  end  event,  representing 
the  achievement  of  the  final  objective ,  at  the  extreme  right  of  the  net- 
work. 

A  path  may  be  defined  as  an  unbroken  chain  of  events,,  with  origin 
at  the  beginning  event  and  terminus  at  the  end  event.   The  network  con- 
sists of  many  such  paths „  overlapping,  paralleling,  and  crossing  each 
other. 

Activities  generally  consume  time,,   The  exact  amount  of  tirae9  t, 
which  will  be  consumed  in  completing  any  one  activity  is  not  generally 
known  in  advance.   If  one  looked  at  the  history  of  a  completed  project,, 
one  could  discover  the  exact  aroount  of  time*  t,  required  to  complete  any 
given  activity,  and  could  label  each  event  with  its  actual  time  of  occur- 
rence.  The  longest  paths  with  respect  to  time,  from  the  beginning  event 
to  the  end  event  could  then  be  found.   This  path  is  called  the  critical 
path.   The  activities  on  this  path  are  called  critical  activities.   The 
length  of  the  critical  path,  the  sura,,  ^£ts  of  the  completion  times  for 
all  activities  on  the  path,  is  the  total  elapsed  time  from  the  beginning 
to  the  end  of  the  project,  and  is  called  the  duration,  D.   It  is  apparent 
that  a  small  change  in  t  for  one  of  the  activities  on  the  critical  path 
causes  the  same  amount  of  change  in  Ds  but  a  small  change  in  t  for  an 
activity  not  on  the  critical  path  would  cause  no  change  in  Do   Hence  it 
is  important  that  project  management  be  able  to  predict  the  critical 
path  in  order  to  predict  Dy  and  in  order  to  allocate  resources s  and  fix 
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schedules,  most  effectively. 

Any  path,  other  than  the  critical  path,,  is  shorter  than  the  critical. 
path.   The  difference  in  Lengths  of  the  critical  path  and  any  other  path 
is  called  the  slack,  K,  for  the  non-critical  path.   An  activity  or  event 
may  be  located  on  several  paths  simultaneously.   The  smallest  slack 
associated  with  any  path  which  includes  a  givers  activity  or  event  is  the 
value  of  slack  associated  with  that  activity /event.   The  activities/ 
events  on  the  critical  path  have  zero  slack.   Any  single  activity  may  be 
delayed  an  amount  less  than  or  equal  to  its  slack  without  affecting  the 
duration  of  the  project. 

The  term,  slack,  may  have  a  slightly  different  meaning  if  some  au- 
thority has  imposed  a  Directed  Completion  Date ,  DCDs,  on  the  end  event  of 
the  project ,   If  so,  we  define  the  Scheduled  Dura t ion ,  D_ ,  as  the  elapsed 
time  from  the  beginning  event  to  the  DCD   Then  slack  for  the  critical 
path  is  (Ds  -  D) .   Slack  for  any  other  path  is  (Ds  -  path  length).   With 
this  definition,  slack  for  any  activity  may  be  either  positive,  negative , 
or  zero,  but  all  activities  on  the  critical  path  will  have  the  same 
slack y  which  will  be  less  than  that  for  any  other  activity  in  the  net- 
work. 

The  PERT  system  is  a  method  by  which  predictions  of  the  critical 
path,  and  slack  for  any  activity*  may  be  made  with  sufficient  accuracy 
to  enable  management  to  operate  more  effectively. 

The  Stochastic  Network  Mode  1 .   in  the  planning  stage,,  before  a  proj- 
ect has  commenced,  the  actual  time,  t,  which  will  be  required  to  complete 
any  given  activity,  is  an  unknown  quantity,   It  is  unlikely  that  this 
time  could  be  predicted  exactly.   We  shall  call  such  uncertain  quantities 
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random  variables,  and  symbolize  thera  with  an  asterisk,,  vs  t*.   The  coi 
spending  unstarred  symbol  will  represent  a  specific  value  of  the  variable 
It  is  possible  to  set  limits,  within  which  t*  is  almost  certain  to  fall . 
It  is  also  possible  to  estimate  the  shape  of  the  probability  density 
function  governing  the  random  variable ,  t*.   The  probability  density 
function  for  t*  is  a  curve  whose  ordinate  in  a  region  is  a  measure  of 
the  likelihood  that  t*  will  fall  in  that  region.   To  be  more  precise, 
the  area  under  the  curve  in  a  region  is  the  probability  that  t*  will 
occur  in  that  region,   in  Figure  3  the  ratio  of  the  shaded  area  to  the 
total  area  under  the  curve  is  the  probability  that  t*  will  fall  between 


t,  and  ty> 


«0 


i 


/. 


x 


tj  t. 

Figure  3 

A  Probability  Density  Function 


The  actual  probability  density  function  for  t*  is  difficult,  if  not 
impossible  to  obtain,  and  will  be  unique  for  any  given  activity  perforate* 
in  a  given  time  period.   Theoretically,  if  the  activity  were  performed 
repeatedly,  under  the  same  conditions 9   with  no  learning  taking  place „ 


then  a  histogram  showing  the  number  of  completions  occurring  in  a  time, 

t  t  At,  plotted  against  t  in  discrete  segments,  would  approximate  the 
probability  density  function,   Since  this  procedure  is  not  possible,  we 
must  be  satisfied  with  the  best  estimate  we  can  formulate,  using  avail- 
able knowledge  of  the  nature  of  the  activity. 

In  order  to  estimate  the  probability  density  function  of  t*  for  a 
given  activity,  the  most  qualified  supervisor  in  charge  of  the  activity 
makes  three  estimates  of  the  completion  time  for  the  activity,   One 
estimate,  $>,  an  optimistic  estimate,  is  the  smallest  amount  of  time  the 
activity  may  be  reasonably  expected  to  consume,.   Another,  M,  is  the 
amount  of  time  the  activity  is  most  likely  to  require,   Finally 8  he 
makes  a  pessimistic  estimate v  P,  which  is  the  greatest  amount  of  time 
the  activity  can  be  expected  to  consume,  barring  completely  unforeseen 
circumstances.   P  and  0   form  the  upper  and  lower  bounds,  respectively, 
of  the  variation  of  t*,  and  M  is  the  mode  of  the  density  function*   The 
probability  density  function  is  assumed  to  have  the  value  zero  at  p  and 
d ,   and  to  reach  a  single  peak  at  M> 

Any  of  several  well  known  mathematical  functions  exhibit  the  prop- 
erties we  have  just  assumed  for  the  distribution  of  t*„   The  Beta 
function, 

f(t)  =  k(t-a)*  (b-t/  (1; 

with  parameters  a.bjcx,  and/  ,  was  chosen  as  the  model  for  the  PERT  sys- 
tem activity  time  because  it  could  be  specified  to  fit  the  constraints 
of  the  estimates.   The  three  estimates  0,  M,  and  P,  specified  three  of 
the  parameters  of  the  Beta  function.   The  remaining  parameter  was  speci- 
fied by  arbitrarily  requiring  the  standard  deviation,  (f  (t*) ,  to  equal 
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Elig  .   One-sixth  of  the  range  is  a  frequently  used  estimator  for  the 
6 

standard  deviation  of  unitnodal  frequency  distributions a  hence  was  con- 
sidered to  be  a  reasonable  assumption.   For  example ,  the  standard  normal 
distribution,  truncated  aft  2.66,  has  its  standard  deviation  equal  to 
1/6  the  range.  [_6j   As  a  result  of  the  above  requirements,  one  of  the 
parameters,  say  <X  ,  becomes  a  function  ©f  the  relative  position  of  M  in 
the  range  between  0  and  P.   Determination  of  this  parameter  requires  the 
solution  of  a  cubic  equation,  (see  Appendix  I).   The  exact  determination 
of  the  mean,  or  expected  value  of  t*,  symbolized  by  te  ,  is  a  tedious 
process,  since  te  is  a  function  of  a,b,o<and  /.  The  following  reasonably 
close  linear  approximation  for  the  mean  was  adopted  as  a  standard; 


E  [t*J  s 


(2) 


The  probability  density  function  thus  assumed  has  a  shape  similar 
to  that  shown  in  Figure  4,   This  distribution  will  be  referred  to  subse 
quently  as  the  PERT  Beta  distribution. 


</ 


M  P 

Figure  4 

A  Typical  PERT  Beta  Density  Function 


A  more  detailed  discussion  of  the  PERT  Beta  distribution  is  included  in 
Appendix  1, 
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For  purposes  of  predicting  the  critical  path,  the  three  estimates 
are  combined  by  means  of  the  weighted  average  formula,  (2),  for  the  mean. 
The  expected  value,  or  mean,  te,  is  used  as  a  deterministic  value  of  t* 
in  the  prediction  of  the  critical  path,   te  lies  between  M  and  the  mid- 
point of  (0,P),  and  is  separated  from  M  by  1/3  the  distance  between  M 
and  the  midpoint.   In  this  manner,  account  is  taken  of  the  usual  tendency 
for  estimators  to  be  optimistic  about  placing  M.   The  fundamental  reason s 
however,  for  using  te  instead  of  M  as  an  estimator  for  t*  lies  in  proba- 
bility theory,  rather  than  an  effort  to  correct  for  this  supposed  bias. 

In  order  to  determine  Slack,  two  new  random  variables,  T*  and  T,*, 
are  defined.   T*  is  associated  with  a  given  event,  and  is  the  elapsed 
time  between  the  beginning  event  and  the  occurrence  of  the  given  event, 
T*  is  the  sum  of  the  random  variables,  t*,  for  the  activities  on  the 
longest  path  to  the  event.   At  this  point  in  the  development  ©f  the 
basic  PERT  system  an  assumption  is  made  that  the  longest  path  to  an 
event  will  be  the  path  whose  mean  length  is  greatest.   The  mean  length  of 
a  path  is  the  sum  of  the  means,  tg,  of  the  activities  on  the  path.   Then9 
by  the  above  assumption,  the  mean,  Tg»  of  T*»  is  the  largest  of  the  mean 
lengths  of  all  paths  to  the  event.  Operating  on  the  basis  of  the  same 
assumption,  the  further  assumption  is  made  that  T*  is  approximately  nor- 
mally distributed.   This  assumption  is  justified  since  we  assume  we  know 
the  longest  path  to  the  event,  and  since  this  path  is  the  sum  of  n  random 
variables,  the  Central  Limit  Theorem  of  probability  theory  asserts  its 
tendency  toward  normality  as  n  increases.   Due  to  the  assumed  independ- 
ence of  the  activities  along  this  path,  the  variance,  d    (T*> ,  of  T*  is 
the  sum  of  the  variances,  (f    (t*),  of  the  activities  on  the  path  to  the 
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event.   The  distribution  of  X*  is  now  completely  specified  by  Tj?  and 
d*2(T*)  in  view  of  the  normality  of  T*. 

We  may  associate  a  random  variable 8  1*^  with  an  activity  as  well  as 
an  event.  We  define  T*  for  an  activity  as  the  elapsed  time  until  the 
activity  is  completed.   Note  that  T*  for  an  activity  is  not  necessarily 
equal  to  T*  for  the  event  following  the  activity.   T*  for  that  event  is 
equal  to  the  maximum  of  the  T*"s  for  the  activities  incident  to  it. 9  that 
is ?  terminating  at  the  event „ 

The  random  variable ,  T^*,  for  an  event  is  defined  as  the  latest 

time  the  event  can  occur  and  not  delay  the  scheduled  completion  of  the 

project.   Tj*  is  calculated  by  subtracting  from  Ds  the  sum  of  the  t*'s 

for  all  the  activities  on  the  longest  path  from  the  event  to  the  end 

event.   Again  the  assumption  is  made  that  this  longest  path  will  be  the 

path  whose  mean  length,  Xc«,»  is  longest.   Then  the  mean  value  ©f  T|*» 

e  x 

written  T\  ,  s  (Ds  -  ^t  ),  where  the  sum  is  the  maximum  of  all  such  sums 
on  paths  between  the  given  event  and  the  end  event.   The  variance  of 
T\*t    d    (T**),  =  d    (t*),  along  the  path  with  greatest  mean  length. 

A  laiddom  variable,  T|*,  may  also  be  associated  with  an  activity. 
Ti*  for  an  activity  is  defined  as  the  latest  elapsed  time  since  the  be- 
ginning event  at  which  the  activity  must  be  completed  in  order  not  to 
delay  the  scheduled  completion  of  the  project*   Notice  that  Tj*  for  any 
activity  is  the  same  as  the  T**  for  the  event  following  the  activity* 

It  is  apparent  that  the  sum  of  the  two  paths  used  to  determine  T* 
and  Tj_*  for  a  given  event/activity  is  the  longest  path  from  the  begin- 
ning to  the  end  which  contains  that  event/activity.   From  our  definition 
of  slack,  we  obtain  the  following  relationship  for  slack,  K*; 
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K*  r  Ds  -  (j*    (Ds  -  Tj*)]  r  T,*  -  T*  (3) 

Then  the  mean  value,  Kg,  of  s  ack  is  the  difference,  T?  -  T  .   Notice 
that  this  expression  for  Kg  depends  upon  the  assumptions  made  in  deter- 
mining T^  and  T„. 

The  Probability  of  Meeting  a   Scheduled  Completion  Date.   We  have 
developed  the  random  variable,  T*,  whose  distribution  we  have  specified, 
by  means  of  simplifying  assumptions,  to  be  normally  distributed,  with 
mean  Tg  and  variance,  <f    (T*) .   The  probability,  Ps,  that  an  event  will 
occur  before  some  scheduled  completion  date  is  simply  P  [jT*  <  T„  I  ,  where 
Tg  is  the  elapsed  time  from  the  beginning  event  to  the  scheduled  comple- 
tion date.   In  particular,  the  probability,  Pj»,  that  the  project  will  be 
completed  before  the  directed  date  is  P  [p*  <  D  ~)  ,   Where  the  duration, 
D*„  is  a  random  variable  equal  to  T*  for  the  end  event.   These  probabili- 
ties can  be  obtained  from  a  table  of  values  for  the  normal  random 
variable. 

The  Predicted  Critical  Path  is  that  path  whose  mean  value  of  slack. 
Kg,  is  smallest.   All  events  or  activities  on  this  path  will  have  the 
same  value  for  Kg,  which  will  be  less  than  the  Kg  for  all  other  events/ 
activities. 

This  completes  the  description  of  the  basic  PERT  system  which  is 
currently  in  general  use.   In  practice,  the  network  representing  the 
project  is  drawn,  events  and  activities  are  assigned  identification  num- 
bers, and  are  assigned  to  specific  organizations  for  responsibility. 
The  supervisors  responsible  for  completion  of  each  activity  then  submit 
the  three  estimates,  0,  M,  and  P.   Several  computer  programs  have  been 
written  which  take  the  activity  number,  preceding  and  succeeding  event 
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numbers,  and  the  three  time  estimates  as  inputs,  then  perform  the  net- 
work calculations.  |_5, 7, 15„20J   The  output  lists  the  activities  in  any 
predetermined  ordering,  giving  Tg,  T^,  Kg,  the  probability,  Pn,  of  meet- 
ing  a  DCD,  c^2(T*),  te,  and  (f~{t*).      The  report  may  list  the  predicted 
critical  path,  and  several  other  piths  in  increasing  order  of  Kp.   Man- 
agement may  take  appropriate  action  by  reallocation  of  resources  if  the 
probability  of  meeting  the  overall  project  DCD  is  unacceptably  small. 
In  order  to  "buy  time"  additional  resources  are  allocated  to,  and  addi- 
tional managerial  attention  devoted  to  those  activities  on  the  predicted 
critical  path,  and  to  those  whose  mean  slack  is  relatively  small.   Re- 
sources may  be  shifted  from  paths  with  greater  Kp  to  paths  with  lesser 
Kg.   In  some  cases  reorganization  of  the  network  is  required v   placing 
more  activities  in  parallel  rather  than  in  series  paths 9  or  eliminating 
activities  whose  products  are  not  entirely  essential  to  the  project. 

The  results  of  original  PERT  computations  may  be  the  basis  for 
establishing  an  overall  project  schedule.   This  schedule  specifies  sched- 
uled completion  times ,  ts ,  for  activities,  and  scheduled  times  ©f  occur= 
rence,  T«,  for  certain  events. 

While  the  project  is  in  progress,  revised  inputs  are  submitted 
periodically ,  showing  percent  of  completion,  and  revised  estimates.   In 
the  case  of  a  partially  completed  activity,  these  estimates  should  be 
considerably  more  reliable  than  those  submitted  before  work  was  commenced. 
These  inputs  are  fed  to  the  computer;,  and  another  network  calculation  is 
run,  as  before,  except  that  completed  activities  are  noted  and  the  actual 
completion  time  is  used  in  the  network  rather  than  an  estimated  time. 
The  results  of  these  runs  may  necessitate  schedule  changes 9  and/or  any 
of  the  actions  previously  mentioned. 
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4„    Critique  of  the  Baste  PERT  System 

The  system  which  we  have  just  described  tends  to  identify  and  call 
attention  to  a  single  path,  the  Predicted  Critical  Path,   This  path, 
though  actually  uncertain,  is  deterministically  established  by  virtue  of 
several  assumptions.   In  most  networks ,  the  activities  on  this  "critical 
path"  will  have  probabilities  of  actually  being  critical  much  less  than 
unity*   Furthermore 9  other  activities s  not  on  this  path,  may  have  proba- 
bilities of  being  critical  larger  than  some  activities  on  the  "critical 
path",.  \JBJ     Hence,  the  system,  by  focussing  attention  on  the  single  path, 
may  deeraphasize  activities  which  are  actually  more  important. 

The  assumptions  which  were  made  in  determining  Tg  and  TL  lead  to 
optimistic  results  in  these  figures,  and  in  Kpo   By  reason  of  the  same 
assumptions,  the  probability  that  an  event  is  completed  before  a  DCD  is 
always  optimistic. 

T*  is  the  random  variable  representing  the  elapsed  time  to  occur- 
rence  of  an  event.  Suppose  there  are  n  different  paths  to  the  event. 
Let  the  random  variables,  P|*,  PV*s  °~"»Pm*»  ^e  £*-®  lengths  of  these  n 
paths,   Now  T*  is  the  maximum  of  these  n  variables,  but  the  maximum  of  a 
set  of  random  variables  is  not  normally  distributed,  even  if  the  indi- 
vidual variables  themselves  are  normal =   in  particular,  T*  does  not  have 
the  same  distribution  as  the  p*  with  the  largest  mean8  nor  does  T*  have 
the  same  mean  as  this  p*.   Furthermore,  the  p.*  are  neither  independent 
nor  identically  distributed,  hence  the  calculation  of  the  distribution 
of  the  maximum  of  this  set  is  so  complicated  as  to  render  it  infeasible 
by  analytic  methods. 

in  general,  the  probability  that  an  event  occurs  before  some 
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specific  scheduled  date  is  the  probability  Chat  none  of  the  possible 
paths  overrun  the  scheduled  date,   By  the  PERT  assumptions,  we  consid- 
ered only  one  path,  namely  that  with  the  greatest  mean  length.   Ob- 
viously, if  there  are  several  parallel  paths  with  nearly  the  same  mean 
length,  then  by  considering  only  the  probability  that  one  of  the  paths 
does  not  overrun,  we  are  being  quite  optimistic  in  our  result,, 

For  example,  suppose  there  are  three  parallel  and  independent  paths, 
each  approximately  equally  likely  to  be  the  critical  path,  and  each  hav- 
ing a  probability  of  being  less  than  Ds  of  about  0.50c   Then  P  v  the 

probability  that  the  project  is  completed  prior  to  the  DCD,  is  the  prob- 

3 
ability  that  none  of  the  path  lengths  exceed  Ds>   Hence,  PD  *  (.5)  s  .125, 

The  PERT  procedure  would  consider  only  the  path  with  the  largest  mean 
length,  and  would  calculate  Pq  s  0.50,   For  networks  with  several  paral- 
lel paths  of  comparable  length,  the  error  would  be  quite  large,  as  in 
the  example,.   In  networks  with  only  one  predominantly  long  path,  the 
error  would  be  negligible. 

In  order  to  obtain  a  more  general  prediction  of  criticaiity  than 
the  single  path  prediction,  we  need  to  calculate  the  probability  that 
each  activity  will  be  on  the  critical  path.   We  shall  call  this  proba- 
bility the  Criticality  Index,  written,  P  ,   Computation  of  Pc  involves 
the  concept  of  determining  the  probability  that  a  particular  pj*  will  be 
the  maximum  of  a  set  of  n  random  variables,  jp, .*,'<>   as  before ,  these  p** 
are  not  independent  nor  identically  distributed,  itence  this  computation 
is  not  feasible  analytically.   It  is  nots  however,  difficult  to  deter- 
mine p£  by  the  Monte  Carlo  method.   This  method  also  readily  yields 
Tg„  T^,  and  Pp  without  resorting  to  the  assumptions  which  were  so 
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troublesome,  but  necessary,  in  the  analytic  approach.  A  detailed  dis- 
cussion of  this  method,  and  suggestions  for  specific  use  of  the  criti- 
caiity  index  will  form  a  substantial  portions  of  this  study* 

-■  PERT  Beta  probability  distribution  was  chosen  as  the  distribu- 
tion for  t*  for  an  activity.   Recall  that  in  specifying  this  distribution 
the  variance  of  the  distribution  was  arbitrarily  chosen  to  be  (F  -  ff>  1  _ 
By  this  assumption,  we  have  specified  the  degree  of  uncertainty  with 
which  the  estimator  makes  his  three  estimates  =   In  other  words  j,  we  are 
denying  the  difference  between  the  case  where  the  estimator  may  predict^ 
with  a  high  degree  of  confidence,  that  t  will  occur  within  a  few  time 
units  of  M  ,  and  the  case  where  t  has  a  higher  likelihood  of  occurring 
near  the  extremes  of  the  range  of  t*«   The  probability  distribution  for 
the  former  case  should  have  a  shape  similar  to  that  showm  in  Figure  59 
while  the  latter  case  would  be  better  represented  by  Figure  6,   It  is 
not  unreasonable  to  assume  that,  under  certain  conditions  supervisors 
will  have  a  basis  for  estimates  with  a  narrow  uncertainty  band  about  the 
mode,  even  though  the  range  of  possible  values  of  t*  might  be  quite 
large.   Under  other  conditions,  the  supervisor  might  be  unable  to  place 
M  with  an  appreciable  degree  of  confidence .   In  the  interest  of  more 
accuracy  in  predicting  the  system  behavior 8  perhaps  a  choice  of  density 
functions  should  be  available  to  the  estimator  to  help  describe  his  con- 
fidence in  his  estimates.   This  would  not  require  a  knowledge  of  proba- 
bility theory  on  the  part  of  the  supervisor.   He  could  simply  be  asked 
to  rate  his  uncertainty  as  to  the  relative  position  of  M  on  a  numerical 
scale,  say  with  three  descriptive  choices.   His  choice  would  lead  to  a 
corresponding  probability  distribution.   Appendices  1^  II ,  III  and  IV 
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describe  nieans  of  providing  this  flexibility  in  choice  of  density  func 
tions  with  the  desired  degree  of  variance. 


Figure  5 

A  Density  Function  with  Small  Variance 
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Figure  6 
A  Density  Function  with  Large  Variance 
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5.    rhe  Monte  Carlo  Technique 

We  have  seen  that  the  completion  time,  t*<,  for  any  activity  in  a 
proposed  project  is  a  random  variable.   As  the  project  is  performed » 
each  activity  is  completed  in  some  actual  time,  t.   Only  after  comple- 
tion of  the  entire  project  can  actual  values  of  t,  slack,  and  elapsed 
time,  T\  for  events  and  activities  be  fixed,  and  the  critical  path  de- 
termined, 

The  Monte  Carlo  method  calculates  average  values  of  these  param- 
eters by  repeatedly  simulating  the  performance  of  the  entire  network, 
and  statistically  averaging  the  results-   The  simulation  of  the  project 
is  accomplished  as  follows: 

A  random  value  of  t  for  each  activity  in  the  network  is  drawn  from 
the  probability  distribution  for  t*  for  the  activity.   These  random 
values  of  t,  called  realizations  of  t*„  are  drawn  in  such  a  way  that  in 
repeated  drawings  from  a  specific  probability  distribution,,  the  number 
of  realizations  falling  within  a  segment  of  fixed  length  on  the  t  axis 
is  proportional  to  the  height  of  the  probability  density  curve  over  that 
segment.   Each  activity  has  its  own  unique  probability  distribution^,  de- 
termined by  the  three  estimates,  hence,  for  one  simulation  of  the  project 
a  realization  of  t*  drawn  for  an  activity  is  a  completion  time  which 
could  have  resulted  from  a  normal  performance  of  the  activity.   One  simu- 
lation of  the  project  consists  of  drawing  a  value  of  t  for  each  activity 
of  the  network.   Calculations  of  the  network  are  then  performed,  using 
these  values  of  t  instead  of  the  te  used  in  the  analytical  method.   T  and 
K  for  each  activity  are  calculated,  and  then  the  critical  path  for  that 
one  simulation  is  determined.   For  each  activity  a  tally  is  incremented 
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If  the  .activity  was  actually  on  the  critical  path.   Another  tally  is  in- 
cremented for  each  scheduled  event  if  the  event  was  accomplished  prior 
to  its  scheduled  completion  date. 

The  above  simulation  Is  repeated  many  times  and  the  values  of  t,  T, 
and  K  for  each  activity  are  added  to  cumulative  totals  after  each  simula- 
tion of  the  project-   After  an  adequate  number  of  replications,  these 
cumulative  totals  and  tallies  are  divided  by  the  number  of  replications, 
yielding  the  average  values,  te  ,  1'   and  Kg  for  each  activity  and  the 
probability  that  the  activity  will  be  on  the  critical  path.   For  all 
scheduled  events,  the  probability  that  the  event  was  completed  prior  to 
the  scheduled  date  is  also  computed.   The  observed  variance  of  these 
statistics  is  also  available  if  desired,  as  well  as  an  approximate  proba- 
bility distribution  for  project  duration. 

The  method  of  sampling  the  probability  distributions  of  the  activity 
times  for  realizations  of  t*  is  discussed  in  detail  in  Appendices  I,,  II, 
III,  and  IV.   Using  the  above  techniques,  it  is  not  difficult  to  provide 
a  choice  of  probability  distributions  for  t*  in  order  to  account  for 
varying  degrees  of  uncertainty  in  estimates,  as  was  discussed  earlier. 

The  computer  codes  designed  to  implement  the  basic  PERT  system  may 
be  modified,,  without  great  difficulty,  to  perform  the  Monte  Carlo  calcu- 
lation.  The  author  recently  modified  the  North  American  Aviation  PERT 
Mark.  Ill  code  in  this  manner.   This  program,  designed  for  the  IBM  7094, 
will  accommodate  a  network  of  1500  activities  without  requiring  external 
memory  units.  Il6j 

The  Monte  Carlo  technique  is  not  subject  t©  the  pitfalls  ©f  inde- 
pendence, assumptions  of- normalcy  due  t©  the  Central  Limit  Theorem,  and 
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assumption  of  a  particular  longest  path,  etc,  which  plague  the  an 
computation.   The  inaccuracies  inherent  in  the  Monte  Carlo  technique  are 
due  to  non- randomness  in  the  random  number  generator „  and  perhaps ,  fail- 
ure to  perform  a  sufficient  number  of  replications  to  permit  the  law  of 
Large  numbers  to  obtain.   Excellent  random  number  gener^t^rs  are  .avail- 
able, providing  resuLts  well  within  the  limit  ©f  accuracy  imposed  by  the 
time  estimates.   Statistical  methods  exist  by  which  to  determine  t 
number  of  replications  required  in  order  to  obtain  the  desired  degree 
confidence  in  the  results.  |_17J 

The  system  is  limited  in  accuracy  only  by  the  human  error  in  estimat- 
ing activity  times ,  and  in  approximating  the  correct  probability  distri- 
bution for  t*. 

The  Monte  Carlo  technique  requires  more  computer  running  time  t 
does  the  analytic  method.   An  efficiently  coded  routine  requires  approxi- 
mately 750  microseconds  to  draw  a  random  number  and  perform  the  calcula- 
tions necessary  to  transform  it  into  a  realization  ©f  t*„   For  details 
of  the  method  employed ,  see  Appendix  1.   The  approximate  increase  in 
PERT  Program  time  required,  if  the  Monte  Carlo  technique  is  used  in  liens 
of  the  standard  analytic  method,,  may. be  estimated  by  the  product? 

(No.  of  activitiesKNo,  of  replications) (750  x  10"6  sec). 

For  example,  a  network  consisting  of  1000  activities  would  require 
12.5  minutes  of  additional  computer  running  time  if  a.  Monte  Carlo  calcu- 
lation of  1000  replications  were  used  instead  of  the  standard  technique. 
The  advantages  inherent  in  this  technique  would  seem  t©  render  this  in- 
crease in  running  time  negligible. 
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6,    he  Extension  to  PERT/COST  LO] 

An  extension  of  the  basic  PERT  system,  PERT/COST,  was  adopted  1  July 
1962  as  the  basic  standard  for  management  time  and  cost  control  systems 
by  DOO  &n<d   NASA.   PERT/COST  extends  the  original  PERT  system  to  provide 
cost  estimation  and  cost  monitoring  features  in  addition  to  the  procedure 
for  prediction  and  control  of  time  performance „   Two  additional  optional 
features,  a  Time-Cost  supplement,  and  a  Resource  Allocation  supplement 
were  also  promulgated, 

PERT/COST  adds  to  the  basic  system  the  following  capabilities; 

a.  Initial  estimation  of  project  cost,  broken  down  by  work  pack- 
ages 9  sub-systems,  etc. 

b.  Provision  of  a  consolidated  estimate  of  the  requirements  for 
manpower  and  other  resources  as  an  aid  to  scheduling  and  pro- 
curement. 

c   Provision  of  in  progress  cost  reports s  showing  current  cost 
status  and  revised  cost  estimates  in  comparison  with  budgeted 
expenditures  and  contract  estimates. 

The  optional  Time-Cost  supplement  provides  a  method  of  estimating 
three  feasible  schedules  for  project  completion,  with  their  estimated 
costs  and  an  evaluation  of  the  technical  risks  associated  with  each 
schedule.   These  schedules  are  designed  to  provide  completion  in  accord- 
ance with  three  criteria;   (1)  the  most  efficient  use  of  time  and  re- 
sources ;  (2)  completion  by  the  DCD;  (3)  earliest  feasible  completion. 

The  Resource  Allocation  supplement  outlines  a  procedure  for  sched- 
uling a  project  to  meet  a  DCD  in  an  optimal  manner  with  respect  to  costs, 
and  utilization  of  manpower, 
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The  concept  of  operation  of  PERT/COST  is  essentially  as  follows; 

^*   Project  breakdown.   The  project  is  broken  down  into  activities 
and  work  packages,  and  graphically  represented  by  a  network.   For  the 
purposes  of  cost  reporting,  estimating,  and  accounting,,  several  activi- 
ties may  be  combined  into  &   "work  package",  when  the  activities  involved 
are  small,  and  detailed  cost  reporting  and  estimating  of  these  activities 
would  involve  unnecessary  expense.   In  general  s,  a  work  package  does  not 
comprise  a.  closed  network.  That  is,  it  may  not  be  integrated  and  repre- 
sented by  a  single  activity  due  to  the  presence  of  events  interior  t©  the 
work  package  which  are  connected  by  activities  to  other  work  packages, 
In  this  case s  the  work  package  must  be  broken  down  into  its  component 
activities  and  events  for  the  network  time  calculations.   The  cost  esti- 
mating and  expenditure  reporting  section  of  the  system  may  utilize  the 
entire  work  package  as  the  smallest  organizational  division. 

Certain  work  packages  may  not  appear  on  the  project  network.   In 
general,  these  are  the  functions  of  management  not  directed  toward  the 
accomplishment  of  specific  activities  and  events.   For  example,,  the 
accounting*  purchasing,  or  management  travel  costs  for  the  project  would 
be  difficult  to  represent  as  activities  in  the  network;  however,  the 
cost  of  these  work  packages  may  be  estimated  and  reported,,  and  treated 
in  the  cost  analysis  portion  of  the  PERT/COST  system. 

B"   Time  calculations.   Time  estimates  are  made  for  each  activity, 
and  the  network  calculations  performed.   The  project  duration  is  com- 
pared with  the  DCD„  and  if  necessary,  any  of  the  following  actions  taken 
in  order  to  reduce  the  overall  project  durations 

(a)  Reorganize  the  network  by  increasing  the  parallelism  of  activi- 
ties, 
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(b)  Commit  additional  resources  to  selected  activities,  or  reallo- 
cate resources  within  the  network. 

(c)  Change  or  delete  activities. 

C.  Manpower  loading.  (9J   A  breakdown  by  skills  of  manpower  required 
for  an  activity/work  package  may  be  included  with  the  time  and  cost  esti- 
mates .   Computer  programs  now  in  existence  assemble  this  information  and 
prepare  a  report  and  display  showing  the  project  manpower  requirements , 

by  skills,  plotted  against  time.  »'7>11J   This  information  provides  project 
management  with  a  forecast  of  manpower  needs ,  and  a  basis  for  reschedule 
ing  activities  in  order  to  best  utilize  available  manpower.   Scheduling 
a  large  project  without  considering  the  overall  use  of  manpower  skills 
leads  to  very  uneconomical  manpower  use.   One  time  period  may  require  a 
manpower  level  for  a  certain  skill  far  in  excess  of  that  available,  while 
the  following  time  period  makes  little  use  of  that  skill   In  performing 
the  necessary  rescheduling,  the  time  constraints  of  the  network  and 
criticalness  of  activities  must  be  considered.   Those  activities  least 
likely  to  be  critical  are  the  obvious  candidates  for  rescheduling  for 
the  purpose  of  leveling  manpower  requirements.   Computer  programs  now  in 
existence,  supplementing  standard  PERT  routines,  perform  this  reschedul- 
ing under  the  constraints  of  maximum  availability  of  manpower  and  with 
consideration  of  slack.  |7,  llj 

D,  Scheduling  the  network.   A  scheduled  duration,  ts,  is  promul- 
gated for  each  activity.   These  times  may  be  less  than,  equal  to,  or 
greater  than  te,  depending  on  the  situation  and  management  policy.   Using 
the  ts  for  each  activity,  the  network  is  computed  and  values  of  Sg  and 
S^,  corresponding  to  TE  and  TL,  determined.   S-  represents  the  earliest 
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date  on  which  an  activity/event  can  be  scheduled  for  completion.   S, 
represents  the  latest  date  on  which  an  activity/event  may  be  scheduled 
for  completion  without  causing  a  schedule  slippage  for  the  project  cons- 
pie  t  ion.  SE  and  S,  are  computed  in  exactly  the  same  manner  as  were  T« 
and  T.  ,  when  t  is  used  as  activity  time  instead  of  t  .   As  a  final  step 
in  scheduling,  an  elapsed  time,  Tg,  is  imposed  at  each  event.   Tg  is 
chosen  to  lie  between  Sp  and  S»  ,  and  may  be  translated  into  a  firm  sched- 
uled  completion  date  for  activities  immediately  preceding  the  event,  and 
a  beginning  date  for  activities  immediately  following  the  event, 

E.  Budgeting  the  project.   After  the  schedule  is  finalized,  cost 
estimates  are  prepared  for  each  work  package ,  and  a  budget  prepared  for 
allocating  funds. 

F.  Monitoring  progress  of  the  project.   As  the  project  progresses, 
reports  are  submitted  from  activities/work  packages ,  indicating  progress, 
new  time  estimates,  funds  committed  to  date,  and  new  estimates  for  cost. 
This  information  forms  inputs  to  the  PERT/COST  computer  progra®,  which 
performs  the  network  calculations,  and  prepares  reports  showing  new 
values  for  Tg,  T»,  KE,  PD,  and  the  new  critical  path,  as  well  as  current 
expenditures  in  comparison  with  estimates  and  budget.   The  new  projected 
cost  curve  is  also  generated  and  compared  with  estimates  and  budget. 
Armed  with  this  information,  management  may  reschedule  reallocates,  re- 
budget,  or  take  whatever  action  as  may  be  appropriate  to  the  situation* 

The  Time -Cost  Option  Supplement.   The  purpose  of  this  supplement  is 
to  standardize  procedures  by  which  a  project  manager  may  prepare  three 
alternative  project  schedule  proposals,  and  evaluate  the  technical  risks 
associated  with  each.   By  technical  risk,  we  mean  the  gamble  of 

26 


performances  time,  and  cost  incurred  by  departing  from  the  best  develop- 
ment technique  in  order  to  meet  the  schedule. 

First  a  plan  is  formulated  designed  to  meet  the  project  BCD,,  as 
imposed  by  the  controlling  agency. 

Nexts  a  Most  Efficient  plan  is  formulated s  in  which  the  best  devel- 
opment techniques  are  used  in  order  to  reduce  technical  risks,  and  most 
efficient  use  is  made  of  manpower  and  resources 9  resulting  in  greater 
duration  for  the  project.   Total  cost  is  generally  lower  for  this  piano 

Finally ,  a  Shortest  Time  plan  is  proposed ,   The  purpose  of  such  a 
plan  is  to  explore  the  feasibility  of  scheduling  the  project  for  a  dura- 
tion shorter  than  that  proposed  by  the  contracting  agency  in  order  to 
benefit  from  any  strategic  benefits  which  might  accrue  from  such  early 
completion.   In  order  to  arrive  at  the  shortest  time  plan„  the  schedule 
is  compressed  by  allocation  of  additional  resources,  paralleling  activi- 
ties which  should  normally  go  in  sequence s  eliminating  activities ,  and 
changing  technical  approaches  to  the  problems  of  design  and  construction. 
All  such  actions  serve  to  increase  the  technical  risks  involved ,  and  most 
of  them  increase  cost. 

The  three  options  are  presented  to  the  contracting  agency  with  esti- 
mates of  cost  and  time,  and  with  a  comparison  and  evaluation  of  the  tech- 
nical  risks  associated  with  each  plan. 

The  PERT/COST  system  is  the  primary  tool  for  obtaining  the  time  and 
cost  estimates  needed  to  prepare  the  three  schedule  proposals. 

The  Resource  Allocation  Supplement ,   This  supplement  provides  a  pro- 
cedure by  which  project  managers  may  schedule  the  project  in  the  optimal 
manner  with  regard  to  the  time-cost  trade-off.   The  concept  of  opera 
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is  as  follows. 

For  any  activity  there  may  exist  several  feasible  schedules,  inv- 
ing  different  time-cost  relationships.   (See  Figure  7.)   Efficient  use 
of  manpower  and  existing  resources  and  machinery  may  result  in  schedule 
A.   Further  time  extension  to  B  might  result  in  higher  costs  due  to  t 
effect  of  fixed  costs.   By  use  of  overtime,  hiring  additional  manpcwer, 
or  purchase  of  additional  machinery  or  space,  time-co'St  combination  C  or 
D  might  be  achieved. 


Cost  A 


t  l  01  ■: 


Figure  7 
Time-Cost  Plans  for  an  Activity 


Schedule  A  would  be  the  most  economical  for  this  activity,  hence 
would  be  the  time  estimated  for  the  activity  for  the  first  network  cal- 
culation,  In  a  similar  manner,  most  economical  estimates  are  obtained, 
with  corresponding  cost  estimates,  from  the  other  activities „   The  net- 
work calculation  is  then  performed,  and  the  project  duration,  D,  compare? 
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with  Ds«   If  D>Dg,  selected  activities  on  the  critical  path  must  be 
shortened  by  resorting  to  a  higher  cost-shorter  time  work  schedule 

For  each  activity  on  the  critical  path,  determine  the  increase  in 
cost,  divided  by  the  decrease  in  time  that  would  occur  in  moving  to  each 
shorter  time-cost  point.   This  is  the  slope  of  the  time  cost  curve  for 
the  corresponding  time  reduction.   The  activity  with  the  lowest  value  of 
slope  is  then  chosen  as  the  activity  to  reschedule. 

If  this  time  reduction  in  the  critical  path  exceeds  the  difference 
in  slack  between  the  critical  path  and  the  smallest  slack  not  on  the 
critical  path,  it  becomes  necessary  to  recompute  the  network  and  deter- 
mine a  new  critical  path.   The  above  process  is  repeated  as  often  as 
necessary  until  D  £.  Dg. 

Slack  paths  are  then  reexamined  to  see  if  any  activities  may  be  ex- 
tended to  lower  cost  points  without  going  critical,   if  fixed  costs  were 
not  a  part  of  the  activity  estimates,  notice  may  be  taken  of  the  fact  thai 
the  most  economical  point  to  operate  is  usually  to  the  left  of  the  mini- 
mum cost  point  with  respect  only  to  direct  costs.   After  determining  the 
most  economical  plan  with  respect  to  direct  costs,  a  test  is  made  to  de- 
termine if  further  time  reduction  below  the  DCD  would  result  in  lower 
overall  costs.   This  is  done  by  adding  on  the  fixed  costs  for  the  project 
and  attempting  further  reduction  in  time.   If  further  time  reduction  re- 
sults in  lower  overall  cost,  then  further  time  reduction  is  performed 
until  total  costs  reach  a  minimum. 

The  schedule  is  then  adjusted  as  necessary  to  reach  an  optimum  lev= 
eling  of  manpower  and  other  resources  within  the  existing  constraints 
and  considering  total  cost. 
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7,    A  Method  for  Resource  Allocation,  using  the  Critical ity  Index 

In  Section  4,  we  pointed  out  the  errors  in  the  concept  of  a  deter- 
ministic prediction  of  the  critical  path.   We  also  showed  the  optimistic 
bias  inherent  in  the  assumption  that  the  longest  path  to  an  event  will 
be  the  path  whose  mean  length  is  greatest.   The  Monte  Carlo  technique 
does  not  depend  on  this  assumption,  hence,  with  this  technique  we  can 
calculate  Pc,  Tg,  and  P  with  accuracy  limited  only  by  the  supervisor's 
estimates  and  our  approximation  of  the  true  probability  distributions  of 
t*  for  the  activities.   We  shall  now  develop  a  system  for  optimum  resource 
allocation  utilizing  the  criticality  index  obtained  from  the  Monte  Carlo 
calculation. 

In  order  to  use  this  method,  we  must  require  that  a  work  package  be 
a  closed  network.   A  closed  network  is  a  network  with  a  single  beginning 
event  and  a  single  end  event.   Activities  exterior  to  the  closed  network 
may  not  be  incident  to,  or  emanate  from,  events  within  the  closed  network 
other  than  the  beginning  or  end  event.   (See  Figure  8.)   In  this  example, 
if  any  activities  not  included  in  the  work  package  should  be  incident  to, 
or  emanate  from  events  B,  C,  or  D,  then  the  network  would  not  be  closed. 
This  same  restriction  is  implicit  in  the  resource  allocation  supplement 
of  the  standard  PERT/COST  system,  although  not  stated  in  the  PERT/COST 
document.   Notice  that  the  closed  network  may  be  integrated  and  repre- 
sented by  a  single  activity  connecting  events  A  and  E.   Some  PERT  com- 
puter programs  now  in  existence  are  capable  of  integrating  such  sub- 
networks for  simplification  of  the  network  for  the  benefit  of  higher 
levels  of  management.  [l5j 
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Figure  8 
A  Closed  Network 


We  next  require  that  the  supervisor  in  charge  of  each  work  package 
submit  as  many  time-cost  plans  as  are  feasible  for  his  work  package  com- 
pletion.  These  plans  are  represented  by  points  on  the  time-cost  curve  in 
the  same  manner  as  previously  discussed  in  connection  with  the  Resource 
allocation  supplement.   No  points  to  the  right  of  the  minimum  cost  point 
are  submitted.   The  time-cost  information  is  represented  as  a  set  of  n 
pairs,  (time, cost),  in  descending  order  of  time  and  corresponding  ascend- 
ing order  of  cost.   For  example: 

14  weeks,  $40,000 

12  weeks,   60,000 

10  weeks,   90,000 
The  time  corresponding  to  the  minimum  cost  point  is  taken  as  the 
most  likely  estimate,  M,  of  activity  completion  time.   The  additional 
time  estimates,  <b   and  P,  are  also  submitted  for  the  minimum  cost  plane 
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These  three  estimates  and  the  corresponding  cost  are  the  basis  of  the 
preliminary  calculation  of  the  network.   From  this  calculation,  the  pro 
ability,  P^,  of  meeting  the  project  DCD,  is  determined.   Recall  that  this 
probability  is  considerably  more  valid  as  a  result  of  the  Monte  Carlo 
technique. 

If  Pq  is  unacceptably  small,  the  network  must  be  compressed  by  re- 
designing the  network  layout  and/or  allocating  more  resources  to  selected 
work  packages  in  order  to  achieve  shorter  completion  time.   If,  after 
whatever  network  redesign  is  considered  feasible,  the  network  must  be 
further  compressed,  the  resource  allocation  phase  of  the  program  is  per- 
formed.  This  phase  consists  of  the  following  logical  steps: 

1.  For  each  time-cost  plan  of  each  activity  having  positive  P^, 
the  following  figure  of  merit  is  computed. 

ETB  =  PC  ^  (4) 

where  ,\t  is  the  reduction  in  time  achieved  by  using  the  T-C  plan  being 
considered,  over  that  used  in  the  previous  network  calculation.  A  C  is 
the  corresponding  increase  in  cost.   ETB  may  be  considered  as  the  expected 
amount  of  time  which  will  be  "bought"  per  unit  cost  by  selecting  the  T-C 
plan  being  considered. 

2.  Form  3-tuples  consisting  of  each  ETB,  its  corresponding  ^  t, 
and  its  activity  identification.   Order  this  list  in  descending  order  of 
ETB, 

3.  Select  from  the  top  of  the  list,  a  sufficient  number  of  reduced 
time  plans  so  that  the  sum  of  these /\t's  does  not  exceed  a  predetermined 
percentage,  p,  of  D.   The  best  value  for  p  may  be  determined  experimen- 
tally.  The  selection  is  made  from  the  top  of  the  list,  but  only  one 

32 


plan,  that  with  the  largest  At,  is  retained  for  any  one  activity. 

4.  For  the  activities  chosen  for  time  reduction,  the  time  coi 
sponding  to  the  new  T-C  plan  is  taken  as  the  new  value  of  M.   The  esti- 
mates for  P  and  0  are  revised  downward  by  multiplying  by  the  ratio, 

M  (new)^   (xhis  procedure  may  be  subject  to  question;  another  method 
M  (old) 

would  be  to  require  P  and  0  estimates  for  each  T-C  plan.) 

5.  With  the  new  time  and  cost  estimates  for  the  accelerated  activi- 
ties, the  network  is  again  calculated,  and  PD  determined.   All  activities 
will  have  a  new  P  .   The  mean  duration,  and  P^  are  again  evaluated.   If 
still  unacceptable,  the  process  may  be  continued  until  either  an  accept- 
able Pn  is  attained,  or  the  budgetary  limit  is  reached.   The  schedule  is 
finalized  on  the  basis  of  the  time  plans  used  in  the  last  network  calcu- 
lation, after  consideration  of  manpower  and  other  resource  leveling. 

By  the  above  method,  resources  are  allocated  where,  probabilistic- 
ally, they  can  be  expected  to  contribute  most  effectively  to  shortening 
project  duration.   The  requirement  for  several  feasible  time-cost  plans 
may  serve  the  additional  purpose  of  forcing  line  supervisory  personnel 
to  consider  the  time-cost  relationship  more  carefully,  allowing  operation 
at  the  most  efficient  point  on  the  time-cost  curve  whenever  feasible. 

An  optimum  value  for  p  can  be  determined  with  a  few  experimental 
runs  of  the  program,   p  should  be  chosen  as  large  as  possible  in  order 
to  decrease  the  number  of  network  calculations  required,  but  small  enough 
so  that  successive  calculations  produce  a  smooth  compression  of  the  net- 
work.  Operation  of  the  program  will  tend  to  decrease  the  differences 
between  the  ETB's  for  the  activities,  which  is  equivalent  to  operating 
the  activities  at  the  same  level  of  expected  marginal  utility.   This  is 
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a  well  known  principle  of  optimization. 

The  inputs  to  the  system,  three  time  estimates  and  a  T-C  function 
specified  by  several  2- tuples,  are  not  complicated,  and  are  compatible 
with  the  concept  of  progress  and  cost  monitoring  envisioned  in  the  PERT/ 
COST  system. 

A  flow  diagram  of  the  logic  used  in  this  resource  allocation  method 
is  shown  in  Appendix  V. 
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8.    A  Military  Application  of  the  PERT  Network  Approach 

One  of  the  most  outstanding  features  of  the  PERT  approach  is  the 
graphical  portrayal  of  the  project  by  means  of  the  network.   By  this  de- 
vice, relationships  between  the  various  tasks  may  be  clearly  visualized*, 
The  insight  thus  gained  by  the  managers  of  the  project,  into  these  com- 
plex relationships,  enables  them  to  function  with  greatly  increased 
effectiveness  in  their  managerial  capacity. 

A  military  operation  is  in  many  ways  similar  to  a  development  proj- 
ect.  The  operation  may  be  divided  into  a  number  of  tasks,  with  many 
interrelationships  existing  between  individual  tasks.   Hence,  the  net- 
work may  serve  the  same  useful  function  in  this  application  as  in  indus- 
trial development. 

The  phases  of  military  planning  known  as  "The  Development  @f  the 
Plan,"  and  "Supervision  of  the  Planned  Action,"  are  most  readily  bene- 
fited by  the  network  approach.   In  the  formalized  military  planning 
sequence,  the  Development  of  the  Plan  phase  has  been  preceded  by  the 
"Estimate  of  the  Situation,"  in  which  the  assigned  mission  has  been 
studied,  possible  outcomes  based  on  enemy  capabilities  and  alternative 
own  courses  of  action  have  been  analyzed  in  a  game  theoretic  matrix, 
and  a  decision  has  been  reached  regarding  a  general  course  of  action  to 
be  pursued. 

The  problem  now  involves  developing  a  complete,  detailed  plan  for 
the  operation,  including  organizing  the  available  forces  and  assigning  to 
each  task  unit  the  appropriate  tasks  which  make  up  the  general  plan,,   At 
this  stage,  detailed  planning  is  done  in  which  training,  acquisition  of 
intelligence,  movement,  communications,  logistics,  and  battle  action  are 
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alt  considered. 

The  plan  begins  with  an  initial  concept,  which  may  be  rer 
by  a  relatively  simple  network,  whose  tasks  are  stated  in  broad  term 
accomplishment.   As  the  planning  becomes  more  detailed,  these  initial 
tasks  may  be  subdivided  into  networks  of  lower-order  tasks.   This  phase 
may  be  done  at  lower  echelons  of  command,  upon  receipt  by  them  of  their 
superior's  directive. 

The  network  will  aid  in  establishing  a  firm  schedule  for  the  op 
tion,  in  much  the  same  manner  as  is  done  in  industrial  projects, 
also  be  a  valuable  aid  for  conducting  briefings  for  subordinates,  since 
it  will  enable  them  to  better  appreciate  the  relationship  of  their  own 
tasks  to  the  whole  operation. 

During  the  operation,  the  network  may  serve  as  a  display,  on  which 
current  progress  is  indicated  as  reported.   The  effects  of  delays  and 
failures  may  be  more  readily  evaluated  and  corrected. 

The  initial  concept  of  a  plan  to  execute  an  amphibious  landing  in 
enemy-held  territory  might  be  represented  by  the  network  shown  in  Figure  9. 

A  network  may  also  be  used  to  lay  out  the  basic  strategy  of  a  pro- 
tracted military  campaign.   In  this  usage,  tasks  would  be  represented  by 
subsidiary  operations,  such  as:   "Seize  Island  B,"  or  "Cut  supply  line 
XY."  The  procedure  is  similar  to  that  employed  in  the  operation  just  dis- 
cussed, except  that  it  takes  place  over  an  extended  period  of  time*   In 
many  cases,  proposed  actions  are  tentative,  or  the  selection  of  alterna- 
tive actions  may  depend  on  future  developments,  such  as  the  success  or 
failure  of  a  previous  endeavor,  or  on  the  subsequent  enemy  reaction,   We 
may  indicate  this  need  for  decision  at  a  particular  point  in  the  network 
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by  showing  a  "decision  box"  in  lieu  of  an  event  at  that  point.   The  deci- 
sion box  has  more  than  one  task  emanating  from  it;  the  decision,  to  be 
made  at  some  Later  date,  determining  which  task(s)  to  be  performed.   The 
network  then  proceeds  to  show  the  planned  action  to  be  taken  for  each 
alternative.   Figure  10  is  an  example  of  a  portion  of  such  a  strategy 
network.  [13 ,  14 J 
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APPENDIX  I 
Monte  Carlo  Techniques  with  the  Beta  Distribution 
The  model  for  activity  completion  time,  t*,  chosen  by  the  originators 
of  PERT,  is  the  Beta  probability  distribution,  characterized  by  its  den- 
sity function:  (2J 

ft*(t)  =  K(t-0)*  (P-t)  ,   0^  t£  P  (5) 

=  0,  elsewhere. 

We  may  reduce  this  to  the  standard  form  of  the  Beta  distribution  by  the 

following  transformation: 

x*  =  t*-  0  (6) 

P  -  0 

The  probability  density  function  of  x*  is: 

fx*(x)  =  Kx  x^(l-x)  ,   0^  t^  1  (7) 

The  variable,  x*,  is  equivalent  to  t*,  where  0=0  and  P  =  1. 

The  two  parameters,  0  and  P,  specify  the  upper  and  lower  limits  of 

the  variation  of  t*.   It  is  now  desired  to  specify  the  remaining  two 

parameters,  od  and  a    ,  so  that  the  mode  of  the  density  function  occurs  at 

M,  and  so  that  the  distribution  has  the  desired  variance. 

Let  ra  »  ^?,  denote  the  mode  of  x*.   By  finding  the  root  of  f'.(x), 
P-0  x 


we  obtain: 


m  =  -T-  (8) 

The  variance  of  x*  is  given  by: 

,  2       (*.+  1)  (f+  1) 

d  (x*>  =  (oc+r+  2)2(oc^r+3)  (9) 

Now  suppose  we  desire  that  the  standard  deviation,  C'Ct*),  be  expressed 

/P-0\  ^>  1 

by/— j— J,  where  d  is  an  arbitrary  constant.   Then  Q   (x*)  =  -r-.   By 
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,2  1 

requiring  that  G  *(x*)  =   -75  ,  we  obtain  the  following  relationship  in  c*6 

Now  c<.  can  be  determined  as  a  function  of  m  and  d  by  the  above  equation, 
With  c<  determined,  if  is  specified  by  solving  equation  (8): 

1    =^ik!E)  (U) 

m 

The  mean  of  x*  is   given  by: 

E[x]    =     ** J  (12) 

Since  >c  must  be  determined  from  the  cubic  equation,  (10),  for  each  value 
of  m,  the  computation  of  E(jx]  is  cumbersome.   In  the  PERT  Beta  models  d 
was  chosen  as  six.   For  this  value  of  d,  Ejx]  is  approximately  linear  in 
m,  and  may  be  approximated  by  the  following  relation: 

E[x]  =  4tn4~  I     .  (13) 

6 

The    transformation  to   t*  yields: 

r  i          0  +  4M  ■+■  P 
Ejt]    =     * Jl !1  .  (14) 

Random  Sampling  from  the  Beta  Distribution.   There  are  several  good 
random  number  generating  routines  which  yield  random  numbers  from  the  uni- 
form (0,1)  distribution.   We  shall  use  the  technique  of  the  Probability 
Integral  Transformation  to  transform  the  random  number,  u,  drawn  from  the 
uniform  (0,1)  distribution,  to  a  corresponding  sample  from  whatever  other 
probability  distribution  we  desire.  |_3J 

This  transformation  depends  on  the  following  theorem: 
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For  any_  random  variable,  v*,  having  the  probability  distribution 
f unc t ion ,  Fv*,  and  the  density  function,  fv*,  define  the  random  variable 
u*  ;  F  ... ,  that  is : 


V" 


u  =  V(V)  =    fv*(t)dt  (15> 

then,  u*  :is  a  uniformly  distributed  random  variable  on  the  interval  (0,1) , 
Hence,  we  may  use  the  inverse  function,  v  r  F  (u) ,  to  transform  the 
sample,  u,  from  the  uniform  (0,1)  distribution  to  the  corresponding  sample, 
v,  from  the  desired  distribution. 
The  distribution  function: 


f     oC  f 

y  (i-y)  dy  (16; 


Fx*(x)  =  KL 

J 


for  the  Beta  variable,  x*,  must  be  calculated  by  numerical  integration. 
In  order  to  utilize  the  procedure  given  above  for  sampling  from  this  dis- 
tribution, we  may  store  the  tabled  distribution  functions  in  the  memory 
of  the  computer  and  use  a  table  look-up  routine  to  enter  the  table  with 
the  random  sample,  u,  and  determine  the  corresponding  value  for  x.   A 
complete  distribution  function  must  be  stored  for  each  of  the  increments 
of  m,  and  for  a  limited  number  of  values  of  d. 

The  input  card  for  an  activity  specifies  the  choice  of  d  which  best 
represents  the  estimator's  uncertainty.   Also  specified  are  the  three 
estimates,  0,  M,  and  P,   With  d  and  M  specified,  the  proper  table  may  be 
selected,   u  is  produced  by  the  random  number  generator.   The  correspond- 
ing value  of  x  is  obtained  from  the  table,  then  transformed  into  a  reali- 
zation of  t*  by  the  following  transformation: 
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t  .  (P-0)x+  0  (17) 

Tables  I  through  V  are  tables  of  distribution  functions  for 
standard  Beta  distribution,  with  d  taking  on  values  of  4  through  8,  re- 
spectively.  A  complete  function  is  given  for  each  of  11  values  of  m. 

As  an  example  of  the  use  of  the  procedure  given  above,  suppose  $,  M„ 
and  P  were  12,  15,  and  20  weeks,  respectively,  and  the  selection,  d  =  6, 
1.  n,ade.  Then  .  .  «£«  ,  .375.   Suppose  we  draw  the  rando*  number, 
.795,  from  the  uniform  random  number  generator.   Entering  Table  ill  for 
m  =  .3,  with  the  argument,  .795,  we  obtain  x  =  .50.   From  the  table  for 
m   s  .4,  linear  interpolation  yields  x  =  .572.   Interpolating  between 
these  two  values  of  x,  we  obtain  x  =  .554  for  m  =  .375.   The  correspond- 
ing realization  of  t*  is  then  8  x -f  12  =  16.4. 

An  efficiently  coded  routine  designed  to  perform  the  above  table 
searches  and  calculations,  requires  approximately  750  microseconds , 

Figures  11,  12,  and  13  are  families  of  curves  of  the  probability 
density  functions  for  x*  corresponding  to  d  s  4,  6,  and  8,  respectively.. 
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DISTRIBUTION    FUNCTION    OF    THE    STANDARD    BETA    FUNCTION    WITH    D    =    4 


A=  .000 

A=  .079 

A=  .174 

A=  .280 

A=  .393 

A=  . 500 

G*  .720 

G=  .714 

G=  .694 

G=  .654 

G=  .589 

G=  . 500 

X 

M=  .000 

M=  .100 

M=  .200 

«=  .300 

M=  .400 

M=  .500 

.05 

.071 

.059 

.0U7 

.035 

.025 

.018 

.  10 

.068 

.062 

.054 

.045 

.037 

.029 

.15 

.233 

.207 

.179 

.  148 

.119 

.093 

.20 

.309 

.280 

.2M7 

.211 

.175 

.  142 

.25 

.381 

.352 

.316 

.276 

.234 

.  195 

.30 

.451 

.421 

.384 

.342 

.296 

.252 

.35 

.516 

.487 

.451 

.407 

.360 

.312 

.40 

.579 

.551 

.515 

.473 

.424 

.  373 

.45 

.637 

.611 

.578 

.536 

.488 

.436 

:SS 

.692 
.743 

.668 

.638 

•  598 

.552 

.500 

.722 

.694 

.658 

.614 

.564 

.60 

.790 

^  .772 

.748 

.715 

.674 

.627 

.65 

.833 

.818 

.797 

.769 

.732 

•  688 

.70 

.872 

.860 

.843 

.819 

.787 

.  748 

.75 

.907 

'   .897 

.884 

.865 

.8  39 

.805 

.80 

.936 

.930 

.920 

.906 

.885 

.858 

.85 

.961 

.957 

.951 

.941 

.927 

.  907 

.90 

.981 

.979 

.975 

.970 

.961 

.949 

.95 

.994 

.994 

.992 

.990 

.987 

.982 

1.00 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

A=  .590 

A= 

.654 

A= 

.694 

A=  .714 

A=  .720 

G=  .393 

G= 

.280 

G= 

.174 

G=  .079 

G=  .000 

X 

M-  .600 

M« 

.700 

M= 

.800 

tf=  .900 

M=1.000 

.05 

.013 

.010 

.008 

.006 

.006 

.10 

.023 

.018 

.015 

.013 

.012 

.  15 

.073 

.059 

.0U9 

.043 

.038 

.20 

.115 

.094 

.080 

.070 

.063 

.25 

.161 

.135 

.1  16 

.  102 

.092 

.30 

.213 

.181 

.157 

.  140 

.126 

.35 

.268 

.231 

.203 

.  181 

.164 

.40 

.326 

.285 

.252 

.227 

.207 

.45 

.386 

.342 

.306 

.277 

.253 

.50 

.448 

.402 

.362 

.331 

.304 

.55 

.512 

.464 

.422 

.388 

.358 

.60 

.576 

.527 

.484 

.4,48 

.415 

.65 

.640 

.593 

.549 

.512 

.477 

.70 

.704 

.658 

.616 

.578 

.541 

.75 

.766 

.724 

.684 

.647 

.610 

.80 

.825 

.789 

.753 

.718 

.681 

.85 

.881 

.852 

.821 

.791 

.756 

.90 

.932 

.911 

.889 

.865 

.834 

.95 

.975 

.965 

.953 

.938 

.916 

1.00 

1.000 

1 

1.000 

1.000 

1.000 

1*000 

Tabl 
A- 

e  I 

f 

(i-y)G  ay 

V 

(x)  . 

f1/ 

(l-y)G  dy 

M   is   the   position   of   the   mode 
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DISTRIBUTION  FUNCTION  OF  THE  STANDARD  BETA  FUNCTION  WITH  0=5 


A=  .000 

A=  .217 

A=  .509 

A=  .870 

A=1.262 

A=1.625 

G=1.826 

G*1.951 

G=2.036 

G=2.  030 

G=1.893 

G=  1.62  5 

X 

M=  .000 

M=  .100 

M*  .2C0 

M«  .300 

M=  .400 

M*  . 500 

.05 

.114 

.078 

.044 

.019 

.007 

.002 

.10 

.104 

.088 

.066 

.041 

.022 

.010 

.15 

.353 

.288 

.210 

.  132 

.072 

.036 

.20 

.455 

.389 

"  .304 

.210 

.129 

.072 

.25 

.546 

.484 

.398 

.296 

.197 

.  120 

.30 

.626 

.570 

.489 

.385 

.275 

.  181 

.35 

.697 

.648 

.575 

.475 

.360 

.251 

.40 

.758 

.718 

.653 

.561 

.447 

.  329 

.45 

.811 

.778 

.724 

.643 

.534 

.413 

.50 

.856 

.8?0 

.767 

.717 

.619 

.500 

.55 

.893 

.874 

.840 

.784 

.698 

.  587 

.60 

.923 

k  .910 

.885 

.841 

.770 

.671 

.65 

.947 

.939 

.921 

.889 

.833 

.749 

.70 

.966 

.961 

.949 

.927 

.887 

.819 

.75 

.980 

.977 

.970 

.956 

.929 

.880 

.80 

.989 

.988 

.985 

.977 

.961 

.928 

.85 

.995 

.995 

.993 

.990 

.982 

.964 

.90 

.998 

.998 

.998 

.997 

.994 

.987 

.95 

1.000 

1.000 

1.000 

1.000 

.999 

.  998 

•  00 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

A=1.894 

A=2.030 

A=2.036 

A=1.951 

A=1.826 

G=1.263 

G=  .870 

G=  .509 

G=  .2  17 

G=  .000 

X 

M=  .600 

M«  .700 

M=  .800 

K=  .900 

M=1.000 

.05 

.001 

.000 

.000 

.000 

.000 

.10 

.005 

.002 

.002 

.001 

.001 

.15 

.018 

.010 

.007 

.005 

.005 

.20 

.039 

.023 

.015 

.012 

.011 

.25 

.071 

.044 

.030 

.023 

.020 

.30 

.113 

.073 

.051 

.039 

.033 

.35 

.166 

.111 

.079 

.061 

.051 

.40 

.230 

.159 

.115 

.090 

.075 

.45 

.302 

.216 

.160 

.  126 

.105 

.50 

.381 

.283 

.213 

.  170 

.141 

.55 

.466 

.357 

.276 

.222 

.185 

.60 

.553 

.439 

.347 

.282 

.236 

.65 

.640 

.525 

.425 

.352 

.296 

.70 

.725 

.615 

.511 

.430 

.365 

.75 

.803 

.704 

.602 

.516 

.444 

.80 

.871 

.790 

.696 

.611 

.532 

.85 

.928 

.868 

.790 

.712 

.632 

.90 

.969 

.934 

.879 

.817 

.742 

.95 

.993 

.981 

.956 

.922 

.865 

1.00 

1.000 

1.000 

1.000 

1.000 

1.000 

Table   II 


ly  a-. 


7)G  dy 
f V   (l-y)°  dy 

■/ft 


M  is  the  position  of  the  mode 
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DISTRIBUTION  FUNCTION  OF  THE  STANOARD  BETA  FUNCTION  WITH  0=6 


A=    .000 

A=    .361 

A=    .891 

A=l.  57U 

A=2.320 

A=3.000 

G=2.872 

G=3.252 

G=3.565 

G=3.  673 

G=3.480 

G=3. 000 

X 

M=    .000 

M=    . 100 

M-    .200 

M=    .300 

M«    .400 

M=    . 500 

.05 

.153 

.090 

.035 

.008 

.001 

.000 

.10 

.131 

.  109 

.071 

.032 

.010 

.  002 

.15 

.449 

.350 

.222 

.  107 

.040 

.012 

.20 

.564 

.472 

.358 

.  195 

.088 

.033 

.25 

.661 

.581 

.454 

.298 

.159 

.071 

.30 

.740 

.675 

.563 

.408 

.248 

.  126 

.35 

.805 

.755 

.662 

.519 

.350 

.  200 

.40 

.857 

.820 

.7*47 

.62U 

.458 

.290 

.U5 

.89e 

.072 

.818 

.718 

.567 

.392 

.50 

.929 

.912 

.874 

.798 

.670 
.761 

.  500 

.55 

.953 

.9M3 

.917 

.863 

.608 

.60 

.970 

N    .964 

.9M9 

.913 

.838 

.710 

.65 

.982 

.979 

.971 

.949 

.898 

.800 

.70 

.990 

.989 

.985 

.973 

.942 

.874 

.75 

.995 

.995 

.993 

.987 

.971 

.929 

.80 

.998 

.998 

.997 

.995 

.988 

.967 

.85 

.999 

.999 

.999 

.999 

.996 

.988 

.90 

1.000 

1.000 

1.000 

1.000 

.999 

.997 

.95 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.00 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

A=3.481 

A=3.672 

A= 

3.565 

A=3.252 

A=2.872 

G=2.321 

G=1.574 

G= 

.891 

G=    .361 

G=    .000 

X 

M=    .600 

M=    .700 

M= 

.800 

M    .900 

M=1.000 

.05 

.000 

.000 

.000 

.000 

.000 

.10 

.001 

.000 

.000 

.000 

.000 

.15 

.004 

.001 

.001 

.001 

.001 

.20 

.012 

.005 

.003 

.002 

.002 

.25 

.029 

.013 

.007 

.005 

.005 

.30 

.058 

.027 

.015 

.011 

.009 

.35 

.102 

.051 

.029 

.021 

.017 

.40 

.162 

.087 

.051 

.036 

.029 

.45 

.239 

.  137 

.083 

.057 

.045 

.50 

.330 

.202 

.126 

.088 

.068 

.55 

.433 

.282 

.182 

.  128 

.099 

.60 

.542 

.376 

.253 

.  180 

.138 

.65 

.650 

.481 

.338 

.245 

.189 

.70 

.752 

.592 

.437 

.325 

.251 

.75 

.841 

.702 

.5M6 

.419 

.328 

.80 

.912 

.805 

.662 

.528 

.421 

.85 

.960 

.893 

.778 

.650 

.533 

.90 

.988 

.956 

.884 

.781 

.665 

.95 

.999 

.992 

.965 

.910 

.820 

1.00 

1.000 

1.000 

1.000 

1.000 

1.000 

Tabl< 

3     III 

T?      .( 

1> 

(1 

-y)G  dy 

<Y  (i-y)G  ay 


M   is   the   position   of   the  mode. 
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DISTRIBUTION  FUNCTION  OF  THE  STANDARD  BETA  FUNCTION  WITH  D  *  7 


.05 
.10 
.15 
.20 
.25 
.30 
.35 
.40 
.45 
.50 
.55 
.60 
.65 
.70 
.75 
.80 
.85 
.90 
.95 
.00 


A=  .000 
G=3.899 
M=  .000 


.189 

.154 

.530 

.651 

.745 

.818 

.874 

.915 

.944 

.965 

.979 

.988 

.994 

.997 

.999 

1.000 

1.000 

1.000 

1.000 

1.000 


A=  .517 
G=4.655 
M=  .100 

.097 

.  128 

.401 

.540 

.658 

.754 

.830 

.886 

.927 

.956 

.975 

*  .987 

.993 

.997 

.999 

1.000 

1.000 

1.000 

1.000 

1.000 


A*1.329 
G=5.316 
M=  .200 

.026 

.071 

.225 

.361 

.497 

.622 

.729 

.815 

.881 

.928 

.959 

.979 

.990 

.996 

.999 

1.000 

1.000 

1.000 

1.000 

1.000 


A=2.  399 
G=5.598 
F=  .300 

.003 
.022 
.083 
.  175 
.292 
.423 
.554 
.675 
.777 
.857 
.915 
.954 
.978 
.991 
.997 
.999 
1.000 
1.000 
1.000 
1.000 


A=3.569 
G=5.353 
M=    .400 

.000 
.004 
.020 
.058 
.124 
.219 
.336 
.465 
.595 
.713 
.812 
.888 
.940 
.972 
.989 
.997 
.999 
1.000 
1.000 
1.000 


A=4.625 
G=4.625 
M=  . 500 

.  000 
.  000 
.004 
.014 
.  039 
.085 
.  157 
.254 
.371 
.500 
.629 
.746 
.843 
.915 
.961 
.986 
.996 
1.000 
1.000 
1.000 


.05 
.10 
.15 
.20 
.25 
.30 
.35 
.40 
.45 
.50 
.55 
.60 
.65 
.70 
.75 
.80 
.85 
.90 
.95 
1.00 


A=5.354 
G=3.569 
M=  .600 

.000 
.000 
.001 
.003 
.01  1 
.028 
.060 
.112 
.188 
.287 
.405 
.535 
.664 
.781 
.876 
.942 
.980 
.996 
1.000 
1.000 


A=5.597 
G=2.399 

M=  .700 

.000 
.000 
.000 
.001 
.003 
.009 
.022 
.046 
.085 
.143 
.223 
.325 
.4M6 
.577 
.708 
.825 
.917 
.974 
.997 
1.000 


A=5.316 
G=1.329 

M=  .800 

.000 
.000 
.000 
.000 
.001 
.004 
.010 
.021 
.041 
.072 
.119 
.185 
.271 
.378 
.503 
.639 
.775 
.894 
.974 
1.000 


A=4.656 
G=  .517 
M=  .900 

.000 
.000 
.000 
.000 
.001 
.003 
.007 
.013 
.025 
.044 
.073 
.  1  14 
.  170 
.246 
.342 
.460 
.599 
.752 
.903 
1.000 


A=3.899 
G-  .000 
M=1.000 

.000 
.000 
.000 
..000 
.001 
.003 
.006 
.01  1 
.020 
.034 
.053 
.082 
.121 
.174 
.244 
.335 
.451 
.597 
.778 
1.000 


Table  IV 


Jo 


f/  (i-y)G  dy 

'A 


M  is  the  position  of  the  mode 
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DISTRIBUTION  FUNCTION  OF  THE  STANDARD  BETA  FUNCTION  WITH  D  =  8 


A*  .000 

A= 

=  .686 

A= 

=1.825 

A=3.  348 

A  = 

=5.010 

A=6. 050 

G= 

=4.916 

G' 

=6. 175 

G  = 

=7.300 

G=7.812 

G  = 

=7.515 

G=6. 050 

X 

M= 

*  .000 

M= 

*  .  100 

M  = 

■  .200 

M=  . 300 

M' 

=  .400 

M=  .500 

.05 

.223 

.101 

.018 

.001 

.000 

.  000 

.10 

.H2 

.145 

.069 

.013 

.001 

.  000 

.15 

.598 

.445 

.223 

.063 

.009 

.001 

.20 

.719 

.598 

.378 

.  154 

.036 

.007 

.25 

.808 

.721 

.534 

.283 

.095 

.024 

.30 

.872 

.815 

.672 

.434 

.192 

.062 

.35 

.918 

.883 

.784 

.  585 

.321 

.  129 

.40 

.949 

.929 

.867 

.718 

.470 

.228 

.45 

.969 

.960 

.924 

.825 

.620 

.356 

:!§ 

.903 

.979 

.960 
.981 

.901 

.751 

.  500 

.991 

.989 

.949 

.854 

.644 

.60 

.995 

1  .995 

.992 

.977 

.924 

.772 

.65 

.998 

.998 

.997 

.991 

.966 

.871 

.70 

.999 

.999 

.999 

.997 

.988 

.938 

.75 

1.000 

1.000 

1.000 

.999 

.997 

.976 

.80 

1.000 

1.000 

1.000 

1.000 

.999 

.993 

.85 

1.000 

1.000 

1.000 

1.000 

^1.000 

.999 

.90 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

.95 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.00 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

A=7.511 

A*7.811 

A=7.030 

A  =  6.  175 

A=4.916 

G-5.007 

G=3.348 

G=1.758 

G=  .686 

G=  .000 

X 

M=  .600 

M-  .700 

M»  .800 

M=  .900 

M*1.000 

.05 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.15 

.000 

.000 

.000 

.000 

.000 

.20 

.001 

.000 

.000 

.000 

.000 

.25 

.003 

.001 

.000 

.000 

.000 

.30 

.012 

.003 

.001 

.001 

.001 

.35 

.034 

.009 

.004 

.002 

.002 

.40 

.076 

.023 

.009 

.005 

.004 

.45 

.146 

.051 

.021 

.011 

.009 

.50 

.249 

.099 

.043 

.021 

.017 

.55 

.380 

.175 

.081 

.040 

.029 

.60 

.530 

.282 

.139 

.071 

.049 

.65 

.679 

.415 

.222 

.  1  17 

.078 

.70 

.808 

.566 

.334 

.  185 

.121 

.75 

.905 

.717 

.470 

.279 

.182 

.80 

.964 

.846 

.624 

.402 

.267 

.85 

.991 

.937 

.777 

.555 

.382 

.90 

.999 

.985 

.904 

.728 

.536 

.95 

1.000 

.999 

.981 

.899 

.738 

.00 

1.000 

1.000 

1.000 

1.000 

1.000 

Table  V  ' 

(/  (1-7)G  dy 
V  U-y)G  dy 


i 


M  is  the  position  of  the  mode 
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Figure  11 
Beta  Density  Functions  .for  d  ~  4 
>s  at  0,  ,2,  and  .5 
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Figure  12 
Beta  Density  Functions  for  d  =  6 
s  at  0S  . 1 ,  . 2 ,  o 3 ,  .4,  and  . 5 
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Figure  13 

Beta  Density  Functions  for  d  =  8 

Modes  at  0,  . 1,  ,2,  „3,  .4,  and  .5 
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APPENDIX  II 
A  Simpler  Method  for  Operating  With  the  Beta  Distribution 
Since  the  Monte  Carlo  method  does  not  require  that  the  variance  of 
t*  be  calculated  and  applied  in  order  to  evaluate  Pp,  we  may  consider  dis- 
pensing with  the  requirement  that  the  variance  of  t*  be  constant  over  the 
range  of  M.   This  requirement  is  not  necessarily  a  natural  one.   The  vari- 
ance of  the  triangular  distribution,  for  example,  defined  naturally  by  0, 
M,  and  P,  is  given  by: 

Cf2(t*)  =  2f?--lf-    (P-M)(M-0)  (18) 

We  see  that  the  variance  for  this  distribution  varies  from  a  maximum  of 

2\iZ.!k]      with  m  at  the  extremes  of  its  range,  to  1  ( ZzE]      with  m  at  the 

1  6  /  216/ 


midpoint.   The  variance  of  the  SinCos  function,  (Appendix  IV),  another 
function  naturally  specified  by  the  three  estimates,  behaves  in  the  same 
general  manner. 

We  may  simplify  the  determination  of  the  Beta  function  parameters, 
avoiding  the  necessity  to  solve  the  cubic  equation,  (10),  by  the  follow- 
ing procedure: 

Fix  the  variance  of  the  symmetric  Beta  function,  (m  =  .5),  by 

Letting 

|  \  2 

_>  9     |  P-0  ! 

a  ^(t*  =    -f-j  d9) 

2 
where  d  is  an  arbitrary  positive  constant  —  \jl2,   or 


Cr2(x* 


=  ™  (20) 

s     d2 

m  =  ,.  j 


For  the  case  where  m  =  .5,  equations  (8),  (9),  and  (20)  yield  the 
following: 
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cL.f.     iLjH 


Define    C  =  $ — lil 

8 


(21) 


(22) 


(23) 


Then,  for  any  value  of  m  in  its  range,  let  cL   and  f   be  determined  as 
follows : 

{    =  C,     0  £  m  ^  .5 
<X  =  C,    .5  £  m  ^  1 
The  remaining  parameter  may  be  determined  by  equation  (8).   The  variance 
of  x*,  as  determined  by  equation  (9),  then  varies  with  ra  as  shown  in 
Figure  11. 


./     -2     -3 


.sr 


,7 


•9    '-o 


Figure  14 
Change  of  Variance  with  Mode  Position 

The  mean,  E[x*]  ,  may  be  calculated  directly  by  equation  (12).   The 
following  close  linear  approximations  for  E  [x*]  ,  for  several  values  of  d 
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the  relative  weights  given  to  0,  M,  and  P  in  the  determination 
expected  value,  hence  give  an  intuitive  feeling  for  the  relationship  be- 
tween d  and  the  amount  of  uncertainty  with  which  m  is  located. 

d  g  [x*3  (24) 

4  20   +-  M  +  2P 

6 


5 

<t> 

+ 

4M 

4- 

P 

6 

0 

+ 

6M 

+ 

P 

8 

8 

By  this  relatively  simple  procedure,  the  parameters  of  the  Beta 

function  for  any  desired  values  of  d  may  be  determined  for  the  purpose  of 

the  Konte  Carlo  calculations.   The  Probability  distribution  functions 

needed  for  the  transformation  from  u*  to  x*,  may  be  extracted  from  a  table 

of  the  Incomplete  Beta  Function.  (jLJ 
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APPENDIX  III 

Use  of  the  Triangular  Probability  Density  in  PERT 
Monte  Carlo  Calculations 

The  triangular  probability  density  function  is  natural Ly  defined  by 
the  three  time  estimates,  0,  M,  and  P.   Since  the  Distribution  function  of 
this  d<  nsity  nay  be  calculated  by  a  simple  formula,  the  transformation 
from  a  sample,  u,  drawn  from  a  uniform  (0,1)  distribution,  to  a  sample 
from  the  desired  triangular  density,  may  be  accomplished  by  a  single, 
simple  calculation.   When  this  procedure  can  be  used,  a  considerable  sav- 
ing in  computer  running  time  is  effected. 

The  technique  is  illustrated  below  for  the  random  variable,  y*,  hav- 
ing the  general  triangular  distribution  with  lower  limit,  0,  upper  limit, 
P ,  and  mode,  M. 
1 


Figure  15 
Triangular  Density  Function 


The  density  function  of  y*  is: 

f   (y)  -    2(y-fl) 
y*      (P-0)(M-0)  ' 

2(P-y) 
"  (P-0)(P-M)  ' 

=   0, 


M  £  y  ^  P 
elsewhere 


(25; 
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The    distribution   function,      F  .,.(>')    s      f  ,v(z)dz,    is 


2 

V(y)  =  T$j%^j  >       6  -  y  -  M 

A     -      9PM   4-     OV\t     _     \t2 


(26) 


p 


-      ±^  +  M   -    2PM -t-  2Py   -   y*    ,   M^  y 
P-0  (P-0)(P-M) 

=      0,      y^0 

=     1,     y^P 

Now   let   u   =   Fy*(y),    and. solve    for   the   inverse    function,    y  =  Fy*    (u), 


and  we   get: 


y  r  0  -h   \/u(P-0)(M-0)    ,  0  e  u  £  *L0  (27 

P-0 


=   P   -     \/(P-M)(P-0)(l-u)    ,      £|  £  u   £  I 

If  we  draw  a  random  number,  u,  from  the  uniform  (0,1)  distribution, 
and  determine  y  by  means  of  the  function,  (27),  then  y*  has  the  desired 
triangular  distribution. 

The  mean,  E(y*J  >  is  given  exactly  by: 

E[y*]  =   0  +  M  +P  (28) 

The  variance,  (f    (y*)»  is  given  by: 

^2(yi)   ,   2^f_   (P-M)(M-6)  (29) 

The  term  on  the  right  varies  from  0,  when  M  is  at  either  extreme,  to 

,P-0\2 


•;.  i-li-)  (  when  M  is  at  the  midpoint. 
V  6 
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APPENDIX  IV 

The  SinCos  Function 

Another  Probability  Function  Suitable  for  PERT 
Monte  Carlo  Calculations 

Another  probability  function  which  is  completely  specified  by  its 
extremes  and  its  mode  is  the  following,  which  we  shall  call  the  SinCos 
function,  for  obvious  reasons. 

The  probability  density  function  is: 


f-*(2) 


ir 


sin 


2(P-0)  **"  2(M-0) 

-H--  cos^-M)   ,   M^z^P 
2(P-0)     2(P-M) 


(30) 


0 


elsewhere 


The  function  has  the  following  shape,  where  M  can  take  any  position  be- 
tween <b   and  P. 


; 


0 


M 


Figure  16 
The  SinCos  Density  Function 
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The  distribution  function,  F  .,.(/. ) 


F,*(z)    =   0, 


fzJt)dt,    is 


J 


-cc 


z* 


s  Md^    fl-cos  Kl^L 
P-0     L  2(M-0) 


z  £  0 

0t   z  <',-  M 


„-  !H*  +  Zd«  sin  IT<'-M1    ,     M  -  2  -  P 
P-0        P-0  2(P-M) 


z  ^  P 


(31) 


Setting  u   =   Fz*(z),   and   solving  for  z,   we  obtain: 

g  ;  0f   2(M-0)   cos'1      1-   (P-0)u  ]      ,      0  ^  u  ^   M-0 
TT  I  (M-0)     I  '    P-0 


x,    ,    2 (P-M)      .    -l 

TT 


P-0  /        M-0'V  M-0  A        *   . 

—  u  -  — r  u     «    — r  *•  u  = 

P-M  ^        P-0j      '      P-0  ~ 


(32) 


The  above  function  transforms  a  random  sample,  u,  from  the  uniform  (0,1) 
distribution  to  the  corresponding  sample,  z,  from  the  SinCos  distribution, 
The  mean,  E  [sTj,  of  the  SinCos  distribution  is: 

E(z*1  =   (IT -2)0+  (4-Tf  )M  +  (1T-2)P  (33) 

Considered  as  a  weighted  average,  we  see  that  the  extremes,  0  and  P,  re- 
ceive more  weight  than  M. 

The  variance,  CT 2(2*)  is  : 

cf2iz1e)    .    A(n--3)(p-0)2  -f  crr2-i6Trf  4O)(p-m)(m-0)  (34) 

In  •»  form  for  numerical  comparison  with  the  PERT  Beta  variance: 

\2 


rf2(*)      =   2-06(^ 


364  ((P-M) (M-0) 


(35) 
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We  see  that  0^(z*)  is  1.88  +  .  18K  times  the  PERT  Beta  variance,  P-6  . 

6 

where  -1^  K  ^  1,  depending  on  the  relative  position  of  M. 

Although,  with  the  SinCos  function,  the  transformation  from  u*  to 
z*  can  be  accomplished  by  a  single  transformation  equation,  computer  run- 
ning time  may  be  as  great,  or  greater,  than  the  table  look-up  routine , 
due  to  the  number  of  operations  required  to  calculate  the  sin"  or  the 
cos   . 
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APPENDIX  V 

FLOW  CHART  FOR  RESOURCE  ALLOCATION  USING  CRITICALITY 

INDEX 


For  each  T-C  plan,  form  the  following  computer  "word". 
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